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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM — 


AN EXPERIMENTAL INVESTIGATION OF REDUCTION IN TRANSONIC 
DRAG RISE AT ZERO LIFT BY THE ADDITION OF VOLUME TO 
THE FUSELAGE OF A WING-BODY-TATL CONFIGURATION 
AND À COMPARISON WITH THEORY 


By George E. Holdaway 
SUMMARY 


Am experimental investigation was made by the free-fall recoverable- 
model technique to assess at zero lift the possibilities of reducing the 
drag-rise coefficients of a wing-body-cruciform-ts1l combination by adding 
volume to the fuselage. The basic features of the test model were an 
unswept &spect-ratio-3.1 thin wing, a fineness-ratio-12.4 fuselage, and 
four 45° sweptback tail surfaces. The tests covered & Mach number range 
of 0.84 to 1.15 with Reynolds numbers of 6,000,000 to 14,000,000, based on 
the wing mean serodynamic chord. 


Considerable reduction in drag-rise coefficient was effected for 
several different modifications by the addition of properly distributed 
volume to the fuselage. In one instance, & reduction in drag coefficient 
was obtained by adding & volume which was almost four times the exposed 
wing volume. The computation method presented in NACA RM A53H17 generally 
predicted the supersonic drag-rise coefficients for each modification 
within 20 percent of the experimental values. As in the above-mentioned 
report, the predictions at a Mach number of one were not accurate. The 
changes in drag-rise coefficients resulting from the modifications were 
generally predicted with better accuracy than the values of drag-rise 
coefficients. 


INTRODUCTION 


During the past year, fuselage indentations of the "area-rule” type 
have successfully reduced the transonic zero-lift dr&g-rise coefficients 
of numerous wing-fuselage combinations. Å summary of the earlier results 
is presented in reference 1. In some cases, where minimum diameters are 
controlled by the engine or other components, fuselage indentation is not 
feasible. Also, for existing aircraft, indentation may be impractical, 
if not impossible. These facts led to the concept of increasing the 
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fuselage volume in proper regions to produce drag reductions comparable 
to those obtained by indentation. An indication that this concept would 
be feasible was obtained independently by an experiment reported in 
reference 1 and by an analysis presented in reference 2. The results of 
reference 1 included a case where the drag-rise coefficient of an air- 
Plane model was significantly reduced by Lengthening the fuselage and by 
adding volume to improve the area distribution of the rearward portion 
of the model; a further reduction was obtained by filling a dip in the 
area distribution for the forward portion of the model. 


The procedure followed in the analytical approach was to use the 
calculation method of reference 2, which is based on the theory of 
reference 3, to determine if reductions of drag-rise coefficients are 
possible with addition of volume to the fuselage, and to determine what 
modifications would indicate sufficient gains to warrant experimental 
investigation. The configuration studied was an aspect-ratio-3.1 unswept 
wing on a fineness-ratio-12.4 fuselage with a cruciform tail. The more 
promising modifications to the fuselage were those designed for minimum 
drag for the configuration at Mach numbers 1.00, 1.05, and 1.14. This 
analysis, presented in reference 2, indicated that addition of volume to 
the fusel&ge would result in substantial reduction in drag-rise coeffi- 
cients, even at supersonic speeds for the M=1.05 and M=1.14 modifications. 


The investigation of this report was undertaken to provide experi- 
mental data for comparison with the predictions presented in reference 2. 
The experimental results would provide additional data for a quantitative 
assessment of the computation method, and would indicate the degree to 
which the reductions in wave-drag coefficients indicated by theory could 
be achieved as measured reductions in drag-rise coefficients. 


The tests were made by the Ames Aeronautical Laboratory at the 
facilities of the Edwards Air Force Base using the free-fall recoverable- 
model technique. The models tested were of large scale resulting in 
Reynolds numbers of 6,000,000 to 14,000,000, based on the wing mean 
aerodynamic chord, for the test Mach number range of M=0.84 to M=1.15. 


SYMBOLS 


CD, zero-lift drag coefficient, er UM IG aha. 


qoy 
t 
Cp,' zero-lift wave-drag coefficient, *oretleal wave drag at zero lift 
Ë qSw 
ADo 
ACD zero-lift drag-rise coefficient, By 


e local chord measured parallel to plane of symmetry 


PU 
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Cy mean aerodynamic chord of the total wing 


AD,  zero-lifií drag rise above subsonic drag level 


H total pressure in the boundary layer 

Ho free-stream total pressure 

M free-stream Mach number 

q free-stream dynamic pressure 

R Reynolds number based on Cy 

S projection of Ss on a plane perpendicular to x axis 


Sg cross-sectional areas formed by cutting the configurations with 
planes perpendicular or oblique to the x axis 


Sw total wing area 


U velocity at the edge of boundary layer 

u velocity in the boundary layer 

x Cartesian coordinate as conventional body axis 

y distance measured normal to the fuselage surface 


boundary-layer displacement thickness 
MODELS 


The dimensions of the unmodified model are given in figure 1, and 
the radii of the fineness-ratio-12.lh fuselage are listed in table I. 
Additional details of the 45° sweptback tail surfaces are given in 
reference 2. The wing used in the investigation was unswept with an 
aspect ratio of 3.1, & taper ratio of 0.39, and a total plan-form area 
of 21.68 square feet. The wing section was elliptical from 0 to 0.5 of 
the local chord and biconvex from 0.5 chord to the trailing edge. The 
maximum wing thickness-to-chord ratio was 3 percent. The wing had no 
twist, dihedral, or incidence, and was of solid aluminum alloy construc- 
tion. The fuselage radii defined in figure 1 are for a minimum-drag body 
of revolution for given volume and length (Sears-Haack body), but behind 
fuselage station 139.4 the theoretical radii and fuselage length were 
extended as dictated by the space required for the recovery mechanisn. 
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The fuselage radii for the three modifications designed to provide 
minimum wave-drag coefficients for Mach numbers 1.00, 1.05, and 1.14 are 
presented in table I. The axial distributions of cross-sectional area 
normal to the longitudinal axis for the basic model and for the three 
modifications are presented in figure 2. The model modified for M=1.00 
is shown in figure 3 to illustrate the comparative size of the model and 
the fact that the changes in radii are quite gradual even though the vol- 
ume added is large. Ail cross sections were maintained circular as in the 
original configuration. Although the general design procedure was pre- 
sented in reference 2, more detailed comments are included in this report 
describing the specific modifications. 


Modification 1, for M=1.00 


Volume was added to the fuselage to alter the normal cross-sectional 
area distribution of the original configuration to that for a Sears-Haack 
body with the same maximum cross-sectional area (fig. h(&)). In this 
case, the values of projected cross-sectional area S are, of course, 
identical with the values of cross-sectional area Sg formed by perpen- 
dicular cutting planes. | 


The type of body shape used for the modification was the same as that 
for the original fuselage (Sears-Haack body; i.e., minimum-drag body of 
revolution for given length and volume), so that the investigation would 
not be affected by an additional variable. The. equation for the body 
radii (fig. 1) differed only in that the maximum radius was increased. An 
additional advantage of the body shape used was that the ends of this type 
of Sears-Haack body are less slender than some other minimum drag shapes 
and would more effectively fair in the bulges in the area-distribution 
curve due to the tail. Modifications were not made behind fuselage sta- 
tion 165 because fuselage indentation would be involved and this was not 
practical because this section contained the recovery mechanism. 


The volume added to the fuselage was 3.63 cubic feet or almost four 
times the exposed wing volume of 0.92 cubic feet. 


Modifications 2 and 3, for M=1.05 and M=1.14 


The design procedure was similar to that used for modification 1, in 
that volume was added to an area-distribution curve to provide & similar 
Sears-Haack shape; however, the procedure differed with respect to the 
type of area-distribution curves used to determine the modification. 


The are&-distribution curves used were average curves based on average 
projected values of Sg obtained with cutting planes tangent to the design 
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Mach cones. The theory upon which this method is based is discussed in 
reference 3. The resultant average area-distribution curves for the 
design Mach numbers are shown in parts (b) and (c) of figure h, with the 
. volume added for esch modification. 


The volume added to the fuselage was approximately three times the 
exposed wing volume for the M=1.05 modification, and twice that volume 
for the 1.14 modification. 


It should be noted that the average area-distribution curves were 
used only in determining the modifications and were not used in computing 
the drag. The individual curves prior to averaging were used to predict 
the wave-drag coefficients. 


INSTRUMENTATION 


Drag measurements were made with two sensitive NACA recording accel- 
erometers which are accurate to 40.0025 g, producing an expected instrument 
accuracy of Cp = 40.0004 at M=1.00 and Cp = 40.0002 at M=1.14. Acceler- 
ometer 1 was located slightly above, and accelerometer 2 slightly below, 
the model center of gravity. 


Pressure measurements were made with a six-cell recording manometer 
which was accurate within £0.05 inch of mercury for pressure readings near 
zero, and was accurate within 2 percent of the full-scale value of 15 
inches of mercury. Mach number was obtained from a calibrated airspeed 
head and was considered to be accurate within M=+0.01. A four-tube pitot- 
pressure rake (fig. 5) was located at fuselage station 100 to measure the 
boundary-layer profile. Tube openings were located about 0.1, 0.3, 0.6, 
and 0.9 inch from the fuselage surface. For two of the tests, base pres- 
sures were determined by manifolding orifices located using an area- 
weighted basis as shown in figure 6. 


All records taken within the model were synchronized by means of a 
1/lO-second chronometric timer. 


TESTS 


The models were released from a carrier airplane at an altitude of 
40,000 feet and allowed to fall freely without propulsion. All surfaces 
were trimmed for zero lift and recovery was initiated at a safe altitude. 
The first two flights were tests of modifications 1 and 2 which were 
designed for M=1.00 and M=1.05, respectively. The third flight was a 
test of the modification for M=1.05 with the tail fairing behind fuselage 
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gtation 190-5/8 cut off to form a flat base. This latter test was made 

to provide data for correlation with possible wind-tunnel tests, and the 
base pressure was measured with four pressure orifices manifolded together 
and located near the center of the base as shown in figure 6(a). 


The last flight was a test of modification 3 for M=1.14. For this 
flight, an effort was made to obtain an indication of the pressure drag of 
the tail falring which also would be of interest in obtaining approximate 
forebody drag for comparison with possible wind-tunnel tests. Seven ori- 
fices were located as shown in figure 6(b) to represent equal portions of 
projected area for a base diameter of 10-1/8 inches. These orifices were 
manifolded together by a large diameter tube (7/8-inch inside diameter). 


The Reynolds number variation with Mach number for the series of tests 
is given in figure T. 


RESULTS 


The experimental results for the three test configurations with tail 
fairing are presented in figure 8. Included in this figure are the theo- 
retical curves of wave-drag coefficients obtained from reference 2. The 
experimental values of subsonic drag coefficients were used to establish 
the datum above which the theoretical wave-drag coefficients were plotted. 
Comparing the experimental and theoretical drag coefficients in this man- 
ner is equivalent to assuming that the level of friction-drag coefficient 
is constant for each modification over the test range of Mach numbers and 
Reynolds numbers. This assumption was considered to be justified because 
& cursory check by available theories indicated that the variation of 
friction-drag coefficient would be of the same order of magnitude as the 
accuracy of the experimental total-drag coefficients. The assumption was 
further justified for the purpose of comparing modifications, since the 
variation of friction-dreg coefficient would be similar for each modifica- 
tion. The tail-fairing drag, presented in figure O(c), was calculated 
using the manifold pressure from the seven pressure orifices which were 
located on an area weighted basis. The experimental results for the modi- 
fication for M=1.05 with the blunt tail are presented in figure 9. Faired 
curves of the experimental data for the three modifications are presented 
in figure 10. Also included in this figure are the experimental data from 
reference 2 for the unmodified configuration. 


The theoretical results for the three modifications and the original 
model, obtained from reference 2, are repeated in figure ll for convenience 
in making comparisons. 


Total-pressure distributions in the boundary layer at fuselage station 
100 are presented in figure 12, for the three modifications, The boundary- 
layer displacement thickness, 5, for each modification was estimated from 


this figure. 


NACA RM A51F22 (Et ce, 7 


DISCUSSION 


A comparison was made between the experimental drag-rise coefficients 
and predicted wave-drag coefficients computed by the method of reference 2. 
The results of this comparison at four Mach numbers are tabulated in 
table II and plotted in figure 13. The supersonic drag-rise coefficients 
from M=1.02 to 1.14 were generally predicted within 20 percent of the 
experimental values. The experimental drag-rise coefficients st these 
supersonic speeds were generally higher than predicted, but this relation- 
ship might vary for configurations other than those tested. At supersonic 
speeds the maximum deviation of theory from experiment was 23.7 percent for 
the unmodified configuration (table IT). The test data for this latter 
case were taken from reference 2 and are not quite as accurate as the test 
data for the three modifications. At a Mach number of one the experimental 
values were always less than the computed values and were poorly predicted 
for all but the M=1.00 modification. 


Of prime interest in this investigation was an evaluation of the 
effectiveness of the modifications to reduce the drag-rise coefficients by 
adding volume. As shown in figure 10, ail the modifications resulted in 
reductions in drag coefficient over the Mach number range of the tests 
despite the fact that they represent additions of volume from two to four 
times the volume of the exposed wing. This result was in accordance with 
the computed results for all cases except for the M=1.00 modification. 
For the M=1.00 modification, & crossover of the drag-coefficient curve 
with that for the unmodified case was expected at M=1.05 (see fig. 11), 
but the experimental data indicated that the crossover would not occur 
until a Mach number of about 1.13, slightly beyond the test range. This 
was traceable to the fact that the drag-coefficient rise for the unmodified 
configuration was larger than predicted and the drag-coefficient rise for 
the M=1.00 modification, above M=1.10, was less than predicted. 


The relative order of drag-coefficient rise for each modification was 
in accordance with the computed results except for the fact that at a Mach 
number of one the M=1.05 modification, even with the cut-off fuselage, had 
a lower drag coefficient than the M=1.00 modification (which should have 
the minimum drag coefficient at this Mach number, as indicated by the com- 
puted results presented in figure 11). This result is attributed to the 
tendency, previously noted in connection with figure 13, for the computed 
values to be least accurate at M=1.00. It would be of interest to study 
this phenomenon by tests of other wing configurations with a fuselage modi- 
fication for M=1.05. 


The quantitative comparison between the computed and experimental 
improvement in drag-rise coefficients effected by the several modifications 
to the original configuration is presented in table III and summarized in 
figure 14. The differences between modifications 2 and 3 are included to 
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illustrate the possibility of experimentally realizing, to a degree, small 
changes in computed benefits. The results show that the computations 
tended to underestimate the benefits due to the modifications by values of 
drag coefficient from 0.001 to 0.002, with few exceptions. Even at & Mach 
number of one, the accuracies (in increments of drag coefficient) with 
which the differences between configurations were estimated tended to be 
better than the accuracies with which the drag coetflclenta of the indi- 
vidual configurations were estimated. | 


Prior to making the tests, it was anticipated that a separation of 
the boundary layer might be caused by the local pressure gradients on the 
body introduced by body shaping; this would introduce drag changes not 
accounted for by the theory. The boundary-layer measurements showed no 
indication of separation even for the M=1.00 modification (which was the 
most severe change) as indicated by the typical boundary-layer velocity 
ratios presented in figure 15. All the profiles obtained indicated that 
the boundary layer was turbulent at fuselage station 100 where the meagure- 
ments were obtained. This is apparent from the agreement between the data 
points and the theoretical curve for turbulent flow. 


SUMMARY OF RESULTS 


This investigation, utilizing tests of free-fall models ai transonic 
speeds to assess at zero lift the possibilities of reducing the drag-rise 
coefficients of an unswept wing-body-tall combination by adding volume to 
the fuselage, has produced the following results: 


1. Considerable reduction in drag-rise coefficient was effected by 
the addition of properly distributed volume to the fuselage. In one 
instance, a reduction in drag coefficient was obtained by adding a volume 
Which was almost four times the exposed wing volume. 


2. The computation method presented in NACA RM A53H17 generally pre- 
dicted the supersonic drag-rise coefficients for each modification within 
20 percent of the experimental values. As in the above-mentioned mn 
the predictions at a Mach number of one were not accurate. 


3. The cbanges in drag-rise coefficient resulting from the modifica- 
tions were generally predicted with better accuracy than the values of 
drag-rise coefficients. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., June 22, 1954 
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Figure 2.- The axial distributions of cross-sectional area normal to 
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Figure 3.- Model modified for M=1.00. The airplane attachment brackets were 
retracted to form a flush surface during free-fall flight. 
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Sears-Haack body, modification 2 
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Figure 4.- Continued. 


17 


18 


120 140 160 180 200 PPO 


for original configuration . 
80 


tt 


E o7 a = 


*ur bg ‘g vere pejos[ouag 


Erie: 


NACA RM A54F22 


Fuselage station, in. 
Figure 4,- Concluded. 


(c) Modification 3, for M=1.14. 
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Figure 6.- Tail configuration and pressure orifice locations for the two drop 
tests in which base pressures were determined. 
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Flgure 7.- Reynolds number variation with Mach number for the tests of the 
modified configurations. 
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theoretical wave-drag coefficients from reference 2 added to the subsonic 
level of the experimental data. 
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(b) Modification 2, for M=1.05. 


Figure 8.- Continued. 
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(c) Modification 3, for M-1.14. 


Figure 8.- Coneluded. 
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Figure 9.- Experimental zero-lift drag coefficients for modification 2, for 
M=1.05, With tail fairing behind station 190-5/8 out off to form a flat 
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Figure 10.- Comparison of the experimental zero-lift drag coefficients for the 
various modifications with the unmodified configuration. 
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Figure 11.- Theoretical computations from reference 2. 
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Figure la. - Total pressure distribution in the boundary layer measured at fuse- 
lage station 100, 
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Figure 12.- Continued. 
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(c) Modification 3, for M=1.14. 
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Figure 13.- Comparison between computed zero-lift wave-drag 
coefficlents and experimental zero-lift drag-rise coeffi- 
clents (see table II). 
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Figure 14.- Comparison between changes in computed zero-lift 
wave -drag coefficients and changes in experimental zero- 
lift drag-rise coefficients due to modifications (see 
table III). 
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(a) Modification 1, for M-1.00. 


Figure 15.- Velocity profile data points in the turbulent 
boundary layer at fuselage station 100. 
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(b) Modification 2, for M=1.05. 
Figure 15.- Continued. 
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(e) Modification 3, for M=1.14. 


Figure 15.- Concluded. 
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